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ABSTRACT 

We present observations and a dynamical analysis of the comet-like main-belt object, (596) Sclieila. 
y-band photometry obtained on UT 2010 December 12 indicates that Scheila's dust cloud has a 
scattering cross-section ~ 1.4 times larger than that of the nucleus, corresponding to a dust mass of 
^ 3 X 10^ kg. V — R color measurements indicate that both the nucleus and dust are redder than 
the Sun, with no significant color differences between the dust cloud's northern and southern plumes. 
We also undertake an ultimately unsuccessful search for CN emission, where we find CN and II2O 
production rates of Qcn < 9 x 10^'^ s~^ and QH2O < 10^^ s^^. Numerical simulations indicate that 
Scheila is dynamically stable for > 100 Myr, suggesting that it is likely native to its current location. 
We also find that it does not belong to a dynamical asteroid family of any significance. We consider 
sublimation-driven scenarios that could produce the appearance of multiple plumes of dust emission, 
but reject them as being physically implausible. Instead, we concur with previous studies that the 
unusual morphology of Scheila's dust cloud is most simply explained by a single oblique impact, 
meaning this object is likely not a main-belt comet, but is instead the second disrupted asteroid after 
P/2010 A2 (LINEAR) to be discovered. 

Subject headings: comets: general — minor planets, asteroids 



1. INTRODUCTION 
Discovered on 1906 February 21, main-be lt asteroid 



(596) Scheila has a diameter of d = 113.34 km (Tedesco et 
[ol. 1^2004) , an orbital period of Porb = 5.01 years, a semi- 
major axis of a = 2.927 AU, an eccentricity of e = 0.165, 
an inclination of i = 14.66°, and a Tisserand parameter 
(with respect to Jupiter) of Tj = 3.209. On 2010 De- 
cember 10.4, observations by S. Larson using the 0.68 m 
Catalina Schmidt telescope showed that Scheila was ex- 
hibiting comet-like activity. Examination of archived 
Catalina data indicated that the object also appeared 
diffuse on December 3, but was point-sou rce-like on Oc - 
tobcr 18, November 2, and November 11 ( Larson|[2010 ). 

If its activity is truly comctary, Scheila would be tlie 
newest member of the small class of objects known as 
main-belt comets (MBCs), which are objects that ex- 
hibit cometary activity {e.g., dust emission) due to the 
sublimation of volatile ices, but occupy stable orbits en- 
tirely confined to the main asteroid belt (Fig.[T]). MBCs 
are unlikely to originate in the outer solar system like 
other comets, a nd are prob ably native to the main belt 
( [Fernandez Tta l. 2002; Hag highipour|2009[ ) . Their activ- 
ity is believed to be triggered by impacts that excavate 



subsurface volatile material (probably wat er ice; Prial 
nik & Rosenberg 2009 Schorghofer '2008"), exposing it 



to direct solar heat ing, drivi ng sublimation and comet- 
like dust emission (Hs ieh fc Jewitt||2006[ ). Noting the 
strong correlation between activity and proximit y to per- 
ihelion for all currently known MBCs, however, Hsieh el\ 



al. (2011b) proposed the alternate possibility that activ- 
ity is in tact primarily dependent on heliocentric distance, 
and only peaks during the post-perihelion portion of each 
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MBCs orbit due to the finite time required for solar ther- 
mal waves to propagate through insulating surface ma- 
terial before reach the subsurface reservoirs of volatile 
material below. Currently, both hypotheses are consis- 
tent with the available evidence, though we note the first 
hypothesis could be ruled out (at least as a universally 
applicable explanation for MBC activity modulation) by 
the determination of a pole orientation for an MBC that 
is inconsistent with seasonal modulation of activity (i.e., 
where peak activity does not occur near a solstice posi- 
tion, as required by the seasonal modulation hypothesis, 
but near an equinox position). Likewise, the second hy- 
pothesis could be ruled out as a universal explanation for 
MBC activity modulation by the discovery of an MBC 
which exhibits peak activity prior to reaching perihelion. 
To date, gas emis sion has never been directly detected 



for a n MBC {e.g., Jewitt et al. 2009 Licandro et al 



2011 1. This present lack of spectroscopic conhrmation of 



outgassing is largely due to the faintncss of MBC activ- 
ity, requiring extremely sensitive observations. Gas de- 
tection efforts are further complicated by the transience 
of MBC activity, requiring that spectroscopic observa- 
tions be secured soon after new MBCs are discovered, 
something that has only been achieved for the more re- 
cent discoveries. Since rapid gas dissipation timescales 
mean that dust emission may continue long after active 
sublimation has actually ceased, even those observations 
may not have been sufficiently prompt to detect gas. 

Despite the lack of direct gas detections, volatile subli- 
mation is inferred to be the driver of activity for all of the 
objects considered to be MBCs. The indirect evidence for 
cometary activity typically includes observations show- 
ing recurrent activity and dust modeling showing contin- 
uous dust emission during each active episode, both con- 
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sistent only with sub limation-driven dust ejection (Hsieh 
a/.||2004 |2011a|b|) 



i'here is however, the case of P/2010 A2 (LINEAR), 
which orbits in the main belt and exhibit ed a comet- 
like dust tail in 2010. [Moreno et al.\ ( |2010| ) initially re- 
ported that numerical dust modehng indicated that the 
object's apparent activity was due to emission that per- 
sisted over several months, implying that it was likely due 
to ice sublimation. However, follow-up modeling based 
on high-resolution Wide Field Camera 3 ima ges from the 



Hubble Space Telescope pewitt et al. 20101 found that 
the apparent activity of 1^/2010 A2 was actually most 
likely a consequence of an inter- asteroid impact, rather 
than the product of sublimation. Observations revealed 
a gap between P/2010 A2's dust tail and its nucleus, 
suggesting that the dust was ejected in a single burst, 
as would be expected from an impact-generated ejecta 
cloud, and was drifting away. T his scenario was corrob - 
orated with further modeling by Snodgrass et al. (2010), 
who employed both ground-based data and data from 
the OSIRIS Narrow Angle Camera aboard the European 
Space Agency's Rosetta spacecraft, and thus were able to 
study the dust tail from two different viewing perspec- 
tives, providing additional constraints on their model. 
Thus, while P/2010 A2 may have appeared comet-like, 
it is likely not a true comet (i.e., an object that exhibits 
sublimation-driven activity) and is better characterized 
as a disrupted asteroid, that is, an object that exhibits 
dust emission as a consequence of an impact (or impacts) , 
and not through the action of sublimating ice. 

All studies of Scheila's 2010 outburst published to date 
( Jewitt et a/1|2011[ iBodewits et aZ."2011'. 'Yang fc Hsieh 



2011| [Ishiguro et at. 2011a b; Moreno et al. 201 lb| con- 



clude that, like for P/2010 A2, Scheila's dust cloud is 
likely be impact-generated, and not a consequence of sub- 
limation. In this manuscript, we describe the results of 
work aimed at better understanding the full nature of 
Scheila and its 2010 outburst that provide additional ev- 
idence supporting this conclusion. 

2. OBSERVATIONS 

We observed Scheila on multiple occasions shortly after 
the discovery of its comet-like activity. Optical imaging 
was performed in photometric conditions on UT 2010 
December 12 using the University of Hawaii (UH) 2.2 m 
telescope, while optical spectroscopy was obtained using 
the 10 m Keck I telescope on UT 2010 December 17. De- 
tails of these observations are listed in Table[lJ Results of 
additional near-infrared observations conducted as part 
of this obser vational campaign are reported in [Yang fc] 
Ssieh| ( |20TI| ). 

Imaging observations employed a Tektronix 
2048x2048 pixel CCD (0'.'219 pixer^) behind Kron- 
Cousins filters. Spectroscopic observations were made 
using Keek's Low Resol ution Imaging Spectrometer 
(LRIS; Oke et al. 1995), a two-channel instrument 



that splits incoming light into red and blue compo- 
nents, and employs dual Tektronix 2048x2048 CCDs 
(0'.'135 pixel"^). We used a I'.'O-wide long-slit mask, 
600/4000 grism, and 560 dichroic, giving a dispersion of 
O.6IA pixel"^ and resolution of ~4A. 

Bias subtraction and flat-fielding were performed for 
UH data, where flat fields were constructed from im- 
ages of the twilight sky. Photometry of Scheila's nu- 



cleus was performed using circular apertures (with back- 
ground statistics measured in nearby regions of blank 
sky to avoid c oma contamination) and calibrated using 
Landolt ( 1992 ) standard stars. 

FUFour spectroscopic observations, to avoid saturation 
from the nucleus while still maximizing the possibility of 
detecting gas emission in the faint dust cloud, we used 
a north-south slit orientation, offset 0'.'6 east of the ob- 
ject's photocenter. North-south dithers (by 30") were 
employed to sample the sky background. We obtained 
two 900s exposures for a total integration time of 1800s. 
Bias subtraction, flat-fielding, and wavelength calibra- 
tion were performed using Low-RedujF software. Object 
identification, extraction, and flux calibration were per- 
formed using the Image Reduction and Analysis Facility 
(IRAF) software package. 

3. RESULTS & ANALYSIS 
3.1. Photometric Analysis 

We find a mean V^-band magnitude for Scheila's nu- 
cleus of mv — 14.22 ± 0.01 mag inside a photometry 
aperture 3'.'0 in radius (Table [2|. This measurement 
represents a photometric excess of Amy = 0.37 mag 
from the quiescent magnitude (my = 14.59 mag) pre- 
dicted by Scheila's lAU ph ase law [Hy = 8.84 ± 0.04, 
G = 0.076 ± 0.06; Warner et aq[20To| , which sig- 
nificantly exceeds the am plitude 61 its measured rota- 
tional lightcurve (0.09 mag; Warner|2006 ) , indicating the 
presence of unresolved near-nucleus corna equivalent to 
^ 40% of the nucleus scattering cross-section inside the 
photometry aperture. 

In a rectangular aperture enclosing the entire dust 
cloud ( "A" in Figure [2]) , we find a total magnitude of 
niv = 13.73 ± 0.02 mag, corresponding to a photometric 
excess of Amy = 0.86 mag. The size of this rectangular 
aperture is chosen to encompass the entire dust cloud as 
observed in U-band before it becomes visibly indistin- 
guishable from the background sky. However, we see in 
Figure |2] that the dust cloud extends somewhat farther 
when observed in i?-band. Measuring the dust cloud's 
i?-band flux using a larger aperture and comparing to 
the fiux measured using the same aperture used for our 
U-band data, we determine that the U-band flux of the 
dust cloud may be underestimated by ~ 10%, and so 
for the purposes of estimating the total mass of Scheila's 
dust cloud, apply this offset to obtain aperture-corrected 
values of my = 13.63 ± 0.02 mag and Amy = 0.96 mag. 
Using the nucleus for calibration, we find a total dust 
scattering cross section of Ad = (1.4 ± 0.3) x 10"* kni^ 
140% of the nucleus cross-section) . Assuming a typi- 
cal bulk d ensity of p = 1640 kg m~ ^ (cf. the Tagish Lake 
meteorite; Hildebrand et a/.|2006 l and an average parti- 
cle radius of = 1 x 10"" m, we can estimate the total 
dust mass, M^, of the cloud using 



-Adpad 



(1) 



Inserting the scattering cross-section we measure for 
Scheila's dust cloud, we thus obtain an approximate dust 
mass of Aid ^ 3 x 10^ kg. This estimate is of course 

^ D eveloped by J. X. Prochaska and available at 
|http: / /www. ucolick.org/ ~xavier/LowRedux / 
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based on numerous assumptions, particularly average 
dust grain size and bulk density, as specified above, and 
is only computed to provide approximate physical con- 
text to the photometric excess that was measured. For 
reference, employing the different g rain size distr i bution 
and density assumptions made by | Jewitt et al. J 2011 1 
(p = 2 000 kg m~^; = 1 x 10^® m) and Bodewits et al. 
|2011[ ) {p = 2500 kg m-3; = 1 x lO""^ m), we would 
instead obtain Md ~ 4 x 10'' kg and Aid 5 x 10^ kg, 
respectively. 



For comparison, Larson (20101 measured Am 



1.24 m ag on 2010 December 3, while [Bodewits et al. 



(2011 1 measured Am y 



0.66 mag on 2010 December 14 

5 and Jewitt et al. (2011) measured Amy = 1.26 mag 
on 2010 December 28 and Amy — 1.00 mag on 2011 Jan- 
uary 5. However, given the different observational and 
instrumental circumstances involved (especially given 
the technical difficulty of observing the extremely bright 
nucleus and comparatively extremely faint dust cloud si- 
multaneously) , with the e xception of the decl ine between 
the two measurements by Jewitt et al. (2011 1, we do not 
consider any dust mass fluctuations implied by compar- 
ing these disparate data sets to be reliable. 

V — R color measurements (Table [2]) of the dust- 
contaminated nucleus indicate that it is redder than the 
Sun, with the entire dust cloud with nucleus flux sub- 
tracted ("DustA" in Tabic [2]) having an effectively iden- 
tical red color, similar to colors measured for other active 



comets and active Centaurs (Bodewits et al. 2011 Je 



witt 



2009). Assuming that the dust contaminating our 



nucleus photometry has a scattering cross-section 40% 
larger than that of the nucleus (estimated above) and the 
same color as the dust cloud as a whole, we find a dust- 
subtracted color for the nucleus of V — R ^ 0.44 mag, 
consistent w ith a D-type spectral classification ( Fornasier 
et al. 2007). To test whether the northern and south- 



ern plumes of Scheila's dust cloud exhibit any compo- 
sitional differences, we measure their colors individually 
("Duste" and "Dustc"), but within estimated uncertain- 
ties, find no significant color differences. 

3.2. Spectroscopic Analysis 
3.2.1. Gas Emission Search 

The most sensitive probe of sublimating gas in a comet 
is CN emission at 3880A. We show the spectral image of 
Scheila from 3700A to 4100A in Figure [sk. In this im- 
age, the horizontal continuum corresponds to reflected 
light from the nucleus. OH sky emission lines are visible 
as vertical bands, and dark Fraunhofer lines in the Solar 
spectrum and prominent Ca H (3933A) and K (3966A) 
absorption lines are also visible. In a two-dimensional 
spectral image, the intensity, /g, of a cometary emission 
lines should be highest near the center of the continuum 
and gradually decrease with increasing distance from the 
continuum, moving in the spatial direction. No spectral 
features near 3880A exhibit such behavior, and we there- 
fore conclude that no gas is detected. 

We also search for CN emission in a one-dimensional 
spectrum extracted from the spectral image using a I'.'O x 
8'.'1 rectangular aperture centered on the continuum. Sky 
background was measured and subtracted using flanking 
regions 10" to 16" from the nucleus. Calibration was 
performed using a nearby flux standard star and a solar 



analog star (Fig. |3b). Shaded regions in Figures [sjD and 
3b indicate where LiN emission is expected, but we again 
nnd no evidence of emission, consistent wit h wo rk by 
Bodewits ef aI](|2011[),|Howefl fc Lovell|(|2011|), andlJehin 
(|2011|. 



et al. 



3.2.2. Gas Production Rate Computations 

To estimate Scheila's CN production rate, we re- 
move the continuum using a scaled solar analog spec- 
trum (Fig. [3b), which should then leave only gas 
emission. Standard errors in three wavelength re- 
gions m the residual spectrum (3760A-3830A; 3830A- 
39OOA, where CN emission is expected; and 3980A- 
4050A; each 70A in width) are 2.7 x lO'^^ 2.0 x 
10~^^, and 1.8 x 10^^^ erg cni^^ s^^ A^^, respectively 
(Fig. [sj:). We choose 2.7 x 10"" erg cm^^ s"^ A'^ 
as a conservative estimate of the uncertainty in the 
CN band. Since 9P/Tempel 1 observations (Meech 



et al] 2011) utilized the same instrumental settings 
used nere, we assume that any CN band in Scheila's 
spectrum will have the same profile as in 9P's spec- 
trum. We therefore find a peak value of any CN band 

of 3 X (2.7 X 10^1^ erg cnr 
1 



A-i), 



10"^'' erg cm"2 s"^ A 

We calculate the integrated CN band flux, /cn, by 
summing the emission flux in the shaded region, obtain- 
ing /cN = 8.8 X 10~^^ erg cm~^ s~^. We then convert 
/cN to the total number of CN molecules, A'cn, using 

LcN = 4^A2/cn (2) 

log A^cN = log £cN + 2.0 log rh - log g{l) (3) 

where A and are in cm and AU, respectively, and 
g{l) is the resonance fluorescence efficiency, which de- 
scribes the number of photons scattered per second per 
radical, in erg molecule"^. During our observations, 
Scheila had a radial velocity of = 2.5 km s~^, for 
which g{l,rh) = 2.7 x 10~^^ erg molecule" "'^ when 



the Swings effect is taken into account ( Schleicher|2010 ) 



Substituting /cn — 8.8 x 10 ° erg cm s ^ , we obtain 



NcN = 5.33 X 10 



,26 



A simple Haser ( 1957) model is used to derive the CN 
production rate, Qcn, trom A^cn, assuming isotropic out- 
gassing, constant radial expansion of the gas coma, and a 
2-step exponential decay process. We use = 1.3 x 10** 
and Ip = 2.2 x 10^ as the effective H aser scale lengths 
at r/i = 1 AU (A'Hearn et al. 1995), and adopt a gas 
velocity of 



1.112 ■ 



-0.4 



km s 



(4) 



( Biver et al.|1997 ). Integrating the computed spatial col- 
umn density model over a rectangle 1800 km x 14480 km 
(the physical size of the extraction aperture in the sky), 
we find QcN < 9 x 10^'^ s~^. Taking average ratios of 
species in previously observed comets (log[(5cN/QoH] = 
-2.5; Qoh/QhsO = 90%) ( jA'Hearn et aL|1995p , we esti- 
mate a water production rate of QH2O < 10^"^ s~^. How- 
ever, given the uncertainties involved in assuming ratios 
of comet species measured at much closer heliocentric 
distances remain unchanged for an object in the main 
asteroid belt, we regard this estimate to be precise, at 
best, to an order of magnitude. 
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3.3. Dynamical Analysis 

To gain a more complete understanding of the circum- 
stances surrounding Scheila's unusual outburst, we also 
consider various aspects of the object's dynamical nature. 
Specifically, we consider its likely origin and whether it 
belongs to an asteroid family. 

To address the first issue, an effort motivated by the 
possibility that cometary objects in the asteroid belt may 
not necessarily originate w here they are currently seen 
(cf. P/2008 Rl (Garradd); pewitt aL| j2009), we per- 
form numerical simulations to assess Scheila's dynamical 
stability. We generate two sets of 100 test particles with 
Gaussian distributions in orbital element space, centered 
on Scheila's JPL-tabulatcd osculating orbital elements, 
where the two sets are characterized by a values equal to 
Ix and 100 X the JPL-tabulated uncertainties for each 
orbital element (Fig. |4^) . We then use the N-body inte- 
gration package. Mercury ( Chambers|1999 ), to integrate 
the orbit of each test particle forward in time for 100 
Myr. 

In three runs using different randomly generated sets of 
test particles, no objects escape from the asteroid belt, 
indicating that Scheila is dynamically stable over this 
time period, and is likely not a recent arrival from else- 
where in the main belt or the outer solar system. We 
note that despite the 100-fold difference in their initial 
dispersions, objects from both the l-a and lOO-fi sets of 
test particles diverge to occupy similar regions of orbital 
element space (Fig. |4|d). This divergence occurs quickly 
(within 10* years) and then remains approximately con 
stant for the 100 Myr test period (Fig 
jects remaining roughly confined to 2. 
2.936, 0.06 < e < 0.28, and 11.5° < i < 18.5°. Objects 
as large as Scheila in the main asteroid belt are often 
simply assumed to be dynamically stable, of course, but 
the example of the Centaur (2060) Chiron shows that 
it is possible for similarly larg e objects to occupy dy- 



sional destruc tion timescales longer than the age of th e 



4h), with all ob- 
9 < a (AU) < 



nami cally unstable orbits {e.g., Nakamura & Yoshikawa 
1993|. As such, given Scheila's unusual comet-like out 



burst, it is useful to have explicit confirmation (provided 
by our simulations) that the object is in fact dynamically 
stable, and is therefore unlikely to be a recently-arrived 
interloper from elsewhere in the solar system. 

Next, given the suggestion that MBCs might be 
preferentially found among members of young aster- 
oid families due to their increased likelihood of hav- 
ing fresh near-surface ice compared to undisrupted as- 
teroids (Hsich 200 9]), we per form a hierarchical cluster- 
ing analysis (Zappala et al. 1990 19941 on proper ele- 
ments for ^135,00 asteroi ds from th e AstDys website 



(http://hamilton.dm.unipi.it/astdys/i Knezevic & Mi- 
|lani||2003i ) to search tor asteroids dynamically related to 
Scheila. We find only six asteroids within a cutoff value 
of 5v' = 120 m s~^, where commonly recognized families 
have tens to thousands of dynamically related members, 
usually within much s maller cutoff values, as shown by 



Nesvorny et al. (2005), despite the fact that their anal- 
ysis was performed when far fewer asteroids were known 
compared to the present day. Such a small number of 
related asteroids indicates that Scheila likely docs not 
belong to a family. 

This conclusion is consistent with Scheila's size: ob- 
jects ^ 100 km in diameter are predicted to have coUi- 



solar system ( [Farinella et a/.||1998} [Bottke et al.\ 2005[ ) 

Thus, any such body observed today is likely to nave a 
primordial origin. The Vesta family, however, shows that 
even extremely large asteroids can be associated with 
families, specifically if those families were pro duced by 



large but non-catastrophic cra tering impacts ( Binzel & 
Xu[ 19931 iThomas et aL||1997D. While such events will 



leave parent asteroids mostly intact, they would still ex- 
cavate large amounts of surface material. If a parent 
asteroid happens to contain ice reservoirs preserved deep 
beneath its surface, by removing much (if not all) of the 
surface material previously insulating that ice against 
sublimation by solar heating, large cratering impacts 
could render that body far more susceptible to subse- 
quent activating impacts by smaller asteroids, such as 
those hypothesized to have triggered activity in other 
MBCs pfaieh fc Jewitt,^2006) . The insignificant number 
of asteroids found to be dynamically related to Scheila 
indicates, however, that not only is Scheila unlikely to 
be a recently-produced fragment of a larger asteroid, it 
is also unlikely to have experienced a large (fragment- 
producing) but non-catastrophic cratering impact in the 
recent past, and so we conclude that its likelihood of 
having recently exposed preserved icy material from its 
deep interior (i.e., several km or more below its surface) 
is low. 

4. DISCUSSION 
4.1. Scheila's Physical Nature 



Yang & Hsieh (2011 1 found that Scheila's near-infrared 
spectrum (from 0.8 to 4.0 /xm) exhibits a consistent red 
slope, has no apparent absorption features, and generally 
resembles spectra of D-type asteroids. Using an intimate 
mixing model incorporating water ice, amorphous car- 
bon, and iron-rich pyroxene, they were able to reproduce 
Scheila's spectrum except for a clearly visible water ab- 
sorption feature that is present in the synthetic spectrum 
but not the observed spectrum. Despite the absence 
of unambiguous evidence of water ice in these observa- 
tions, given the signal-to- noise ratio of the data, |Yang fc| 
Hsieh (2011 ) concluded that the presence of water ice on 
Scheila s surface could not be excluded to a level of a few 
percent. They further noted that the similarity between 
Scheila and D-type asteroids, which in turn have been 
noted for their similarities t o classical cometary nuclei 
from the outer solar system ( Fitzsimmons et al.| |l994), 
suggests that Scheila c ould contain pres erved ice deep 
within its interior {e.g., Jones et aL||1990 ). 

For comparison, two other MBCs, 133P/Elst-Pizarro 
and 176P /LIN EAR, have been identified as B-type or F - 
type asteroids (Bagnulo et al. 2010[ Licandro et aL|2011 ), 
where both types are subgroups or C-type asteroids. In- 
terestingly, these spectral classifications mean that, un- 
like Scheila, they are more physically simil ar to nearby 
m ain-be lt asteroids than to classical comets ( Licandro et 
a/.||2011 1 , even though we believe the active episodes tor 



these objects were actually sublimation-driven in nature, 
while Scheila's was not. The nuclei of MBCs 238P/Read, 
P/Garradd, and P/2010 R2 (La Sagra), as well as that 
of P/2010 A2, have yet to be characterized spectr oscopi 



cally due to the small sizes of 23 8P and P/2010 A2 ( |Hsieh 
et aL|20l"Tb Jewitt et aL|2010 ) and significant cometary 
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activity at the time of all currently published observa- 
tions of P/Garradd and P/La Sagra ( Jewitt et a/.||2009 



Moreno et a/.]r2011a| . 

As discussed above (Sectionfll, direct spectroscopic de- 
tections^ gas emission in MEKJs have b een elusive (e.g. 



Jewitt et 2009 ■ ,Licandro et aZ. |2011 ). Successful un- 
ambiguous detections ot ice have been likewise difficult 
to obtain. No spectroscopic evidence of exp osed water 
ice has bee n found on any of the MBCs (e. g 
et al. 2011 ). However, Hsieh & Jewitt (2006 1 hypothesize 



that cometary activity tor these objects is being driven 
by small areas of exposed subsurfac e ice (hundre ds of 
square mete rs on km-sized bodies; cf. |Hsieh et a/.||2004[ 
Hsieh 2009). If that approximate ratio of active to in- 
active surface material is correct, we would actually ex- 
pect unresolved disk-integrated reflectance spectroscopy 
to be of limited use without sensitivity levels to one part 
in ~ lO'* or better. 

While this sensitivity level is beyond the reach of cur- 
rent Earth-bound facilities for the km-scale (a nd smaller) 
MBCs at the dist ance of the ast e roid b elt, [Rivkin fc 
Emery (2010) and Campins et al. (2010) have reported 
water ice detections (correspondiiig to rnuch larger sur- 
face coverage than predicted for the MBCs) for the 100- 
km-scale asteroid (24) Themis, which belongs to the same 
Themis asteroid family that also contains 133P and 176P 
(Hsieh & Jewitt"2006^. A similar absorption feature has 



also beerT detected on outer belt asteroid (65) Cybele 
(Licandro et al. 2011). The attribution of the absorp- 
tion feature observ ed by these group s to water ice has 
been challenged by Beck et al. (2011) though, who sug- 
gest that the absorption feature can be equall y well ex- 
plained by th e non -volatile mineral goethite. Rivkin fc 
Emery (2010) and Campins et al. (2010) acknowledge 
the thermal instability ot water ice that they claim to 
detect and propose various scenarios how such surface 
ice could be maintained. However, none of these scenar- 
ios has yet been confirmed to actually plausibly account 
for a widespread, long-lived surface layer of water ice as 
implied by their observations. We further note that no 
outgassing or dust emissio n has ever been ob served for 
(24) Themis (also noted by [Rivkin fc Emery|r2010) , and 
as such, the connection between ice on mam- belt aster- 
oids and the activity of MBCs remains unsubstantiated. 

In summary, while the hypothesis that activity in 
MBCs is sublimation-driven is supported by indirect ev- 
idence such as numerical modeling re sults and observa- 



tions showin g recurrent activity {e.g., Hsieh et al. 2009 
2010 2011b I, it remains unsubstantiated by the imam 



biguous detection of either ice or sublimation products. 
As such, if the non-detections of ice or gas emission via 
spectroscopy currently do not preclude the designation 
of an object as an MBC, the same non- detec tions of ice 
and gas emission for Scheila (Section 3.2 ,Howell fc Lovell 
20Tn IJehin et a/.| [20lT| [Bodewits k al.\[2i}ll\ [Va^IgT 
Hsieh 201 1[ ) necessarily must also be considered inad- 
equate criteria for concluding that the object is a dis- 
rupted asteroid, at least on their own. We must also 
consider other evidence. 

4.2. Morphology Analysis and Outburst Scenarios 

The most obvious aspect of Scheila's dust cloud that 
should bear clues as to its origin is its unusual morphol- 
ogy. Unlike many other comets which typically exhibit 



a single tail, often pointed in the antisolar direction, 
and a coma, Scheila exhibits two distinct curved dust 
plumes extending to the North and the South (Figure [2]). 
Deeper and higher resolution im agery also shows a faint 
westward-pointing dust "spike" ( Jewitt et al.^2QV\\ . The 
two main dust plumes are not easily explained by the 
standard cometary dust ejection paradigm where grains 
are ejected isotropically and subsequently f ollow well- 



Rousselot defined syndyn e and synchrone curves {cf. Finson & 



Probstein 1968), and therefore suggest the action ot a 
more unusual ejection m echanism 



Ishiguro et al. (2011b) present numerical modeling re- 
sults that show that the cloud's morphology can be ac- 
counted for by a hollow cone of dust ( as expected from 
an impact; e.g., Richardson et aL|2007 ) ejected Sunward 
and then turned back by radiation pressure. The hollow- 
ness of the cone is suggested by the dust cloud's apparent 
limb brightening (i.e., Scheila's northern and southern 
plumes), consistent with the greater optical depth ex- 
pected of such a structure along its edges. A solid cone 
of ejected dust (as expected from a sublimation-driven 
jet) similarly pushed back would instead exhibit central 
brightening due to greater optical depth in the jet's core, 
and not appear to exhibit multiple plumes. The dis- 
parate strengths of the northern and southern plumes 
may indicate an oblique angle of incidence for the im- 
pact, with m ore dust expected down range of the inbound 
impactor (cf. [Ish iguro et aL||2011b| 



While the scenario mbdeied by Ishiguro et al. (2011b I 
plausibly accounts for the appearance of multiple dust 
plumes with a single impact, if Scheila's dust cloud is 
produced by a sublimation-driven process, any scenario 
that similarly accounts for multiple observed dust plumes 
likely requires multiple active sites. If these active sites 
are coUisionally excavated, multiple impacts would have 
to have occurred on timescales shorter than their deple- 
tion timescales. This scenario is unlikely, however, given 
the typically low rate of impacts on any single body in 
the asteroid belt (e.g., Farinella & Davis 1992) and the 
expected short depleti on timescale s tor surface volatiles 
on main-belt asteroids (Hsieh 2009') . Multiple active sites 
could be possible if near-surface ice is abundant, trigger- 
ing sublimation at multiple points via thermal stresses. 
This scenario is contradicted though by Scheila's history 
of observed inactivity until 2010, and the fact that it first 
exhibited a comet-like dust cloud aXv^ 240° , well before 
perihelion, when the surface was receiving only ^ 60% 
of the solar flux that it would at perihelion. If Scheila 
were particularly icy, we would expect to have observed 
past activity and at times of maximum solar heating. We 
therefore find any scenario invoking multiple sublimation 
sites to be implausible. 

One possible scenario which could explain the observed 
dust morphology via a sublimation-driven phenomenon 
would be one in which directed dust emission (i.e., a 
jet) is being produced from a single active site near the 
rotational pole of the object, which in turn is directed 
approximately Sunward. Such a situation would mimic 
the effect of a hollow cone of material (i.e., the rotating 
jet) being ejected Sunward, where the cone would then 
be expected to be pushed back by radiat ion pressure as 
in the i mpact-driven model described by [Ishiguro et al.\ 
(2011b). Asymmetry in the density of material forming 
the cone could then be caused by diurnal effects where 
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sublimation, and therefore dust emission, increases in in- 
tensity when the active site is at its closest to the subsolar 
point and receives maximum Solar heating, and decreases 
when the active site is farthest from the subsolar point 
and it receives somewhat diminished Solar heating, as- 
suming the rotational pole is not directed exactly at the 
Sun. 



Given Scheila's rotation period of 15.848 hr (Warner 
Ishiguro et aL|[2011a l, this scenario would require 



2006 



Schella to have thermal inertia low enough for the tem- 
perature of the active site to vary enough during each 
rotation such that the dust production rate of the jet 
visibly varies as the active site rotates into and out of 
direct Solar heating. This scenario would also require 
a somewhat fortuitous rotation pole orientation where 
the angle between the pole and the Sun-object vector 
is large enough to create significant diurnal temperature 
variations over the course of a single rotation, but small 
enough that the projected appearance of a Sunward- 
pointing hollow cone of ejected material that is then 
pushed back by radiation pressure is not lost. The ther- 
mal and numerical dust modeling required to test this 
scenario is beyond the scope of this paper. We note, 
however, that the obh que impact scenario described by 
Ishiguro et al. (2011b) is far simpler than one involving a 
rotating, diurnaliy-varying cometary jet oriented exactly 
in such a way that it mimics the asymmetric Sunward- 
directed hollow cone of ejected material that otherwise 
arises naturally in the impact scenario. As such, in the 
absence of compelling evidence that sublimation must be 
taking place, we favor that impact scenario as the most 
plausible explanation for Scheila's unusual dust emission, 
and therefore consider it to be the second disrupted as- 
teroid to be discovered, after P/2010 A2. 

4.3. MBCs vs. Disrupted Asteroids 

New and upcoming all-sky survey facilities like the 
Panoramic Survey T elescope And Rapid Response Sys- 



2010) and the Larg 



Jones et al. 2009 



tem (Pan-STARRS; [Kaiser et al. 
Synoptic Survey Telescope (LSS' 

as well as long-running existing surveys like LINEAl 
Spacewatch, and the Catalina Sky Survey, are expected 
to discover many more comet-like objects in the future, 
some of which may be found to orbit in the main aster- 
oid belt. A natural question to pose, then, is how can we 
distinguish with any degree of certainty which of these 
objects are MBCs, where comet-like activity is due to 
sublimation, and which are disrupted asteroids, where 
comet-like activity is due to a recent impact? 

Unfortunately, a definitive a nswer to this q uestion does 
not yet exist. For this reason, Jewitt et al. 
taken the approach of designating ail comet 



201 1[ ) have 
ike objects 

with main-belt orbits as MBCs, based on their appear- 
ance alone (i.e., using the observational definition of a 
comet as any solar system object exhibiting mass loss in 
the form of extended surface brightness features such as a 
coma or tail). Using this definition, they classify Scheila 
as an "impact-activated MBC" , agreeing with this work 
in terms of physical conclusions, if not in terminology. 
We prefer, however, to retain the physical meaning of a 
"comet" as an active body whose gas or dust emission is 
due to the sublimation of volatile ices. As such, we choose 
to define MBCs as only those comet-like objects with 
main-belt orbits whose activity is sublimation-driven (ac- 



cording to our best jud gement), as originally specified in 
Hsieh fc Je\^ ( |2006[ ), and objects like P/2010 A2 and 
Scheila, where dust emission is the direct consequence of 
an impact, as disrupted asteroids. 

Of course, imposing the requirement that activity must 
be sublimation-driven for an object to be considered an 
MBC implies that we know the underlying cause of ap- 
parent activity for every comet-like main-belt object, 
which is not necessarily always the case (or indeed ever 
the case, at least not until the first successful gas detec- 
tion is made). 

From a practical standpoint, the development of an un- 
ambiguous method of distinguishing sublimation-driven 
dust emission and impact-driven dust emission will re- 
quire the discovery and thorough individual investiga- 
tions of many more examples of both types of objects 
before a set of reliable criteria can be developed for each, 
meaning that significant additional work is required be- 
fore such a classification system can be achieved. 

In the meantime, we can draw the following prelimi- 
nary conclusions and observations regarding the differ- 
ences between MBCs and disrupted asteroids, based on 
the currently extremely limited populations of each: 

1. Long-lived activity is not a reliable indicator of a 
dust emission event's source. All five currently 
recognized MBCs exhibited activity persisting over 
weeks or months ([Hsieh et al. 2004 , 2009, 2011a][je^ 
witt et a/.[[2009| IMoreno et al., 2011a), but so dJd 



the twocurrently r ecogn i zed disrupted asteroids 



(Jewitt et al. 
Long-lived ac 



2010 2011 



Snodgrass et al. 2010) 



ivity in the five recognized MBCs is 



regarded as an indication of continuous emission 
over an extended period of time, a characteristic 
of cometary emission. However, similar long-lived 
"activity" (>1 month), or rather the persistent 
presence of a dust cloud, in P/2010 A2 and Scheila 
are explained as a consequence of the impact events 
on each object ejecting large 1 /Ltm) particles, 
which sir nply take longer to dissipate than smalle r 
particles ( Jewitt et aL|2010"l Bodewits et a/.|20TI ). 



Observations showing steadily increasing quantities 
of dus t during long-lived active episodes {e.g., for 



238P; Hsieh et al. 2011b) could be better indica- 



tions of sublimation-driven activity since they can 
not simply be explained by persistent large dust 
grains. Ongoing dust particle production is re- 
quired. However, while this can result from the 
actual release of new dust from the nucleus, an in- 
crease in visible scattering surface area could also 
arise via fragmentation of existing dust particles, as 
likely occ urred during Corn et 17P/Holmes's 2007 
outburst ([Hsieh et a/.[2010|. 



In the case of Scheila, [Jewitt et al.\ (f2011[) argued 

the 



that observations of rapid fading of the dust cloud 
(approximately one month after its discovery) was 
an indication of impact-driven emission. While the 
short time between the appearance of Scheila's dust 
cloud (in December 2010) and this rapid fading 
(observed in January 2011) means that this inter- 
pretation is likely correct, we note that decreasing 
activity (with no earlier observations of increas- 
ing activity) was also observed for 133P in 2002 
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dHsieh et aL | | 2004[) and P/Garradd in 2008 (|Je 



witt et al. 2D09p . Thus, while increasing activity 



may suggest cometary emission, and decreasing ac- 
tivity may suggest impact-driven emission, neither 
should be considered conclusive proof of either sce- 
nario without other supporting evidence. 

2. Dust clouds with unusual morphologies may in- 
dicate unusual fo rmat ion circumstances. As dis- 
cussed in Sect ion 4.2 (and referen ces within) for 



Scheila and by Jewitt et al. (^2010 1 and Snodgrass 



et al. (20101 for F/2U10 A2, the unusual structures 



of their respective dust clouds appear to be indica- 
tions of non-sublimation-driven origins. In the case 
of P/2010 A2, an observed gap between the nucleus 
and the dust tail strongly suggested that a dust 
cloud had been ejected in a single impulsive event 
and was slowly drifting away from the nucleus, con- 
sistent with being an ejecta cloud produced by an 
impact, and inconsistent with cometary dust emis- 
sion. The implications of Scheila's pair of dust 
plumes are considered at length in Section |4.2[ 
While such interpretations can be used to draw 
preliminary conclusions, we caution against rely- 
ing solely on morphology-based determinations of 
a dust emission event's source. Multi-tail struc- 
tures have been observed for comets believed to 
exhibit subli mation-driven activity (e.g., 238P and 
P/La Sagra; |Hsieh et aZ. |2009| [Moreno et al.\2011a ) 
and the abrupt termination of sublimation-driven 
dust emission that could cause an observable gap 
between a nucleus and a dust tail, while proba- 
bly uncommon, could be possible given certain to- 
pographical circumstances {e.g., shadowing of iso- 
lated active sites). Detailed numerical dust model- 
ing must be employed to ensure that any unusual 
morphological structures cannot be plausibly ex- 
plained by both sublimation-driven and impact- 
driven dust emission. 

Recurrent activity, separated by intervening periods 
of inactivity, is extremely difficult to explain as the 
consequence of impacts. Upon the initial discov- 
ery of a comet-like main-belt asteroid, conclusions 



the result of numerical dust modeling (|Boehnhardt 


et al. 19961 


Hsieh et aL||2004 2009 2011a| Moreno 


et a/.|20lla 


. The underconstrained nature of these 



are usually never unique, and simply reflect the 
most plausible ejection scenarios as judged by the 
authors of those studies. 

In contrast, observations of recurrent activity for a 
single object le ave far less room for ambiguity (c/. 
133P and 238P; [Hsieh et aL|2004[ [20Tol |201 lb[ ) . In 
those cases, repeated active episodes occurred on 
timescales of ~5-6 years, which correspond closely 
to the orbital period of each object, and are far 
shorter than the expected timescales of repeated 
random collisions on individual main- belt objects 



by impactor s of appreciable size {e.g., Farinella & 



Davis 1992). One might envision a scenario m 
which a main-belt object routinely encounters an 
overdensity of impactors at certain parts of its or- 



bit, perhaps due to an intersecting cometary debris 
stream, for example. The contrived nature of such 
a scenario (in which debris streams dense enough to 
reliably impact 133P and 238P on each orbital pas- 
sage only intersect these specific objects and do not 
cause comparable events on other asteroids with 
similar orbits) argues strongly against its plausi- 
bility though. On the other hand, repeated ac- 
tive episodes are routinely observed for comets for 
which activity is driven by sublimation, and can 
be naturally explained for the MBCs as the con- 
sequence of seasonal variations in solar heating of 
isolated active sites, or perhaps simply from the in- 
crease in solar heating experienced by ea ch object 
near perihelion ( [Hsieh et a/^[2004[ [201l"b[ ). 

As discussed above, the sizes of the known populations 
of both MBCs and disrupted asteroids, from which these 
generalizations are drawn, are extremely small, and so 
the typical caveats associated with small-number statis- 
tics certainly apply. However, for the above reasons, we 
suggest that repeated activity is the least ambiguous and 
most reliably obtainable indicator available at the cur- 
rent time that comet-like activity is sublimation-driven 
for a particular object. While only two of the seven 
currently known comet-like main-belt objects {i.e., both 
MBCs and disrupted asteroids) have been observed to ex- 
hibit recurrent activity to date, we note that the remain- 
ing five objects were discovered as comet-like bodies re- 
cently enough that they have not actually yet completed 
full orbits since their respective discoveries. As such, 
continued monitoring of all of these objects to search for 
recurrent activity will be important for validating their 
identification as MBCs or disrupted asteroids. 

5. SUMMARY 

We have conducted a photometric, spectroscopic, and 
dynamical study of comet-like main-belt asteroid (596) 
Scheila, and report the following findings: 

• Scheila's dust cloud comprises a F-band photomet- 
ric excess of A = 0.96 mag over the expected qui- 
escent brightness of the nucleus, giving it a scat- 
tering cross-section ~1.4 times larger than that 
of the nucleus, and an estimated total mass of 
Afn ~ 3 X 10"^ kg. Separate V - R color mea- 
surements of each of the two curved plumes of the 
dust cloud show that they are both redder than the 
Sun, with no signficant color differences between 
the plumes themselves. 

• We find no evidence of CN emission in Scheila's 
dust cloud down to the sensitivity level of our 
spectroscopic observations, and place upper lim- 
its to CN and H2O production rates of Qqn = 
8.91 X 10^3 s-i and Qh^o 10^^ s'K 

• An analysis of Scheila and its dynamical vicinity 
using numerical simulations indicate that it is dy- 
namically stable for > 100 Myr, suggesting that it 
is likely native to its current location. We also find 
that it does not belong to a dynamical asteroid 
family of any significance, meaning that it is un- 
likely to be the product of a recent fragmentation 
or cratering event that could have exposed deeply 
buried ice. 
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Considering the available evidence obtained by us and 
others, we conclude that Scheila is most likely a disrupted 
asteroid, and not an MBC (where dust emission is the 
result of the sublimation of ice, as in other comets), in 
agreement with previous work. This conclusion is largely 
founded on consideration of the unusual morphology of 
Scheila's dust cloud, coupled with the simplicity of the 
impact scenario required to reproduce it and the inabil- 
ity of any physically plausible sublimation-driven scenar- 
ios to do likewise. While the currently known examples 
of MBCs and disrupted asteroids suggest that morpho- 
logical analyses can provide strong indications as to the 
nature of comct-likc activity in main-belt objects discov- 
ered in the future, wc caution against relying on such 
evidence on its own. We suggest instead that recurrence 
(or absence of recurrence) of comet-like activity is the 
best indicator available to date of whether an object is 
likely to be a MBC or a disrupted asteroid, underscoring 



the importance of long-term monitoring observations of 
these objects. 
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TABLE 1 

Observation Log 



UT Date 


Telescope 


Obs.^' 






Filters 




1/° 




AS 






2010 Dec 12 


UH 2.2m 


Im 


81 


1770 


V 


0"8 


239.6 


3.107 


2.529 


16.4 


-1.0 


2010 Dec 12 


UH 2.2m 


Im 


97 


2230 


R 


0"8 


239.6 


3.107 


2.529 


16.4 


-1.0 


2010 Dec 17 


Keck I 


Sp 


2 


1800 




I'.'O 


240.5 


3.100 


2.462 


15.6 


-1.5 



^ Type of observation (Im: imaging; Sp: spectroscopy) 
^ Number of images 

Total effective exposure time in seconds 

FWHM seeing in arcsec 

True anomaly in arcsec 
' Heliocentric distance in AU 
s Geocentric distance in AU 
^ Solar phase angle in degrees 

' Orbit plane angle (between the observer and object orbit plane as seen from the object) in degrees 



TABLE 2 
Nucleus and Dust Photometry 





my 


rriR 


V -R 


Nucleus^ 


14.22±0.01 


13.77±0.01 


0.45±0.02 


TotaU'' 


13.73±0.02 


13.27±0.02 


0.45±0.03 


DustA'^ 


14.84±0.02 


14.36±0.02 


0.48±0.03 


DustB"* 


22.79±0.03 


22.78±0.03 


0.51±0.04 


Dustc" 


23.43±0.04 


23.00±0.04 


0.43±0.06 



Photometry of nucleus (including unresolved near- 
nucleus coma) inside a 3'.'0 (radius) circular photom- 
etry aperture 

^ Photometry of nucleus and dust inside aperture "A" 
as marked in Figure [2] 

Nucleus-subtractodphotometry of dust inside aper- 
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Fig. 1. — Plots of eccentricity (top) and inclination (bottom) versus semimajor axis showing the distributions in orbital element 
space of main-belt asteroids (black dots), main-belt comets (red circles), and likely disrupted asteroids P/2010 A2 and Scheila 
(blue circles). Also marked with dotted lines are the semimajor axes of Mars (aMars) and Jupiter (ojup), the semimajor axis of 
the 2:1 mean-motion resonance with Jupiter, and the loci of Mars-crossing orbits and Jupiter-crossing orbits. 
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Fig. 2. — Composite ii-band image of Scheila constructed from data obtained using the UH 2.2 m telescope on 2010 December 
12 and comprising data equivalent to 2230 s in total exposure time. Rectangular photometry apertures used to measure dust 
fluxes are marked, where aperture "A" is 53" x 62" in size, "B" is 17"5 x ll'.'O, and "C" is 9"2 x 15'.'9. 
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Fig. 3. — (a) Long slit spectral image of Scheila, taken on UT 2010 December 17 using Keck I. The wavelength range (3830A < A < 
3905A) where CN emission is expected is marked. Sky emission lines are visible along the spatial direction (the vertical axis) and the solar 
Ca H and K absorption lines are marked, (b) Sky background-subtracted and flux-calibrated spectrum of Scheila (solid black line). The 
scaled spectrum of the solar analog, G2V star HD 91163, is marked as a blue dashed line. The shaded region indicates where CN emission 
is expected, (c) Spectrum of Scheila with the underlying solar-type continuum subtracted. Red error bars show the la uncertainties in the 
three wavelength regions discussed in the text. 
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Fig. 4. — Plots of semimajor axis versus eccentricity (top) and inclination (bottom) for (a) the initial orbital elements and (b) final orbital 
elements after a 100 Myr dynamical integration of 100 l-cr Scheila test particles (red dots) and 100 100-cr test particles (blue dots), and (c) 
the orbital elements of the same 100 Icr test particles shown at the start of our simulation (overlapping black dots at the center of each 
panel) and 100 l-cr test particles after 20 Myr (purple), 40 Myr (blue), 60 Myr (green), 80 Myr (orange), and 100 Myr (red). 



